Contacts Dynamics across Jamming in a granular glass experiment 



C. Coulais,^ R. P. Behringer,^ and O. Dauchot^ 

^SPHYNX/SPEC, CEA-Saclay, URA 2464 CNRS, 91 191 Gif-sur-Yvette, France 
^Department of Physics and Center for Nonlinear and Complex Systems, 
Duke University, Durham, North Carolina 27708-0305, USA 
^EC2M, ESPCI-ParisTech, UMR Culliver 7083 CNRS, 75005 Paris, France 

We study experimentally the vicinity of the Jamming transition by investigating the statics and 
the dynamics of the contact network of an horizontally shaken bi-disperse packing of photo-elastic 
discs. Compressing the packing very slowly, while maintaining a mechanical excitation, we produce 
a granular glass, namely a frozen structure of vibrating grains. Nevertheless, we observe in this 
glass phase a remarkable dynamics of the contact network^ which exhibits strong dynamical hetero- 
geneities. Such heterogeneities are maximum at a packing fraction 0*, distinct and smaller than 
the jamming packing fraction which is indicated by the abrupt variation of the average number 
of contact per particle. We further demonstrate that the "dynamical transition" taking place at 0* 
originates in the spatial correlation of the coordination number. Our findings are discussed in the 
light of recent numerical and theoretical studies of thermal soft spheres. 
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At large packing fraction, disordered packings of par- 
ticles with repulsive contact interactions jam into a rigid 
state. For frictionless and a-thermal particles the jam- 
ming transition coincides with the onset of isostaticity 
and a number of geometrical and mechanical quanti- 
ties exhibit clear scaling laws with the distance to jam- 
ming [1]. One prominent signature of jamming is the 
singular behavior of the average number of contacts per 
particle z — zj ex {(j) — where zj = 2(i, d being 

the space dimension ^ |3| . The distribution of the gaps 
between particles displays a delta function at zero and a 
square root decay for increasing gaps, which is key to the 
singular behavior of the average contact number [4j-t7]. 

Although the average coordination number singular- 
ity is the hallmark of jamming at zero temperature, its 
behavior is less clear at finite temperature. Both ex- 
perimentally [8 and numerically, f6\ ISUTT] it was ob- 
served that the first peak of the partial pair correlation 
functions has a finite maximum at a packing fraction 
^j(T) > ^j(O). This maximum has been interpreted 
as a vestige of the divergence seen in the pair correlation 
function at point J, the T = jamming transition. The 
vicinity of point J has been explored in a mean-field-like 
replica description of thermal soft and hard spheres [12 . 
This description recovers all the observed scalings in tem- 
perature and packing fraction but the square root singu- 
larity of the pair correlation function when T = 0+ and 
(j) = . This discrepancy, together with the onset of 
a diverging length in the vibrational properties of the 
jammed state, [13 suggest that larger scale correlations 
must be taken into account, and calls for a better char- 
acterization of the vicinity of point J. 

In the present letter, we focus on the dynamics of the 
contact network, a natural quantity of interest as soon 
as T > 0, which to our knowledge has not been consid- 
ered so far. To this aim, we experimentally investigate 



a vibrated two-dimensional bi-disperse packing of grains 
close to its jamming transition. We verify that the struc- 
ture, in terms of particle neighborhoods, is completely 
frozen on the experimental time scales, hence justifying 
the appellation of a granular glass, and we demonstrate 
the emergence of heterogeneities in the dynamics of the 
contact network. These heterogeneities are maximum at 
a packing fraction (/)* smaller than the jamming pack- 
ing fraction ^j, which is signaled by an abrupt variation 
of the average number of contacts per particle. We fur- 
ther demonstrate that the " dynamical transition" taking 
place at 0* originates in the spatial correlation of the 
coordination number. 

Our experimental set up (Fig. [l] top) is adapted 
from [14", TF, and allows us to shake photoelastic discs 
between two cross-polarizers, thereby accessing the con- 
tact and force network [3l [16]. A bi-disperse mixture 
of 4964 and 3216 polyurethane (PSM-4) discs, of respec- 
tivley 4mm andSmm diameter, lies on a flat transparent 
surface, which oscillates with an amplitude of 1 cm at a 
frequency of 10 Hz. The grains are confined in a static 
rectangular cell, for which the position of one horizontal 
boundary is tuned using a micro-metric piston attached 
to a force sensor. An LED back-light and a circular po- 
larizing sheet are inserted in the vibrating surface, so that 
the grains are lighted by transmission. A monochromatic 
high resolution CCD camera takes two types of images 
in phase with the vibration. Every odd period, a circu- 
lar analyzer is inserted by a rotating wheel in the cam- 
era field to visualize the photoelastic pattern. From the 
direct-light images acquired every even period, we ob- 
tain the grain positions with a resolution of 0.5% of the 
mean grain diameter. Standard tracking and tessellation 
techniques are then used to obtain the dynamics and the 
structure of the packings. From the cross-polarized im- 
ages, we estimate the normal force between two neighbors 
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FIG. 1: Experimental setup, (color online) Top: The vi- 
brating cell. Bottom: The photo-elastic pattern at contact. 
White circles are the grains, red lines the link between neigh- 
bors, and dark lines the Delaunay triangulation. 

by integrating the square spatial gradient of the light in- 
tensity over the area defined by the two Delaunay trian- 
gles sharing a common edge (see Fig. [l] bottom). Hence- 
foth, all lengths are expressed in units of the small grains 
diameter, and time is expressed in units of vibration cy- 
cles. Note that as compared to the brass discs, used in 
[Ml [15] , the photo elastic discs are much softer (Young 
modulus of 0.5GPa, as compared to lOOGPa). 

In order to ensure the highest and most reproducible 
jammed packings, the packing fraction, which we control 
within a relative resolution of 5 10~^ is increased by steps 
of 5(j) = 3.10"^, using exponentially increasing time steps 
(Fig. 11 (top left)). All images are then acquired during 
stepwise decompression; following a decompression, we 
fix the packing fraction, and then carefully check that 
both the pressure - defined as the force measured at the 
wall normalized by the weight of the grains - and the 
dynamics are stationary. Fig. [2] (top right) displays the 
pressure with respect to the packing fraction. Ptot (re- 
spectively Pstat) is the pressure measured when the vi- 
bration is switched on (respectively off). Pstat is thus 
the static pressure sustained by the packing in the ab- 
sence of vibration, whereas Pdy n = Ptot — Pstat can 
be interpreted as the additional dynamic pressure, in- 
duced by the vibration. At high packing fraction, the 
pressure is dominated by its static part, whereas at low 
packing fraction, it is mostly dynamic. A crossover oc- 
curs at some intermediate packing fraction, which in [T4] 
was shown to coincide with a maximum in the hetero- 
geneity of the dynamics, as probed at scales of the or- 
der of 10~^ grain diameters. As compared to the pre- 




FIG. 2: Obtaining a granular glass (color online) Top 
Left : Logarithmic increase followed by a stepwise decrease 
of the packing fraction, and pressure measured at the wall. 
Top Right : Average pressure vs packing fraction: (Q) • 
Ptot, (□): Pstat, (A): Pdyn as defined in the text. Left: 
Superposition of the Laguerre cells computed at times t = 1 
and t = 5000 for the loosest packing (0 = 0.8031). Right: 
Average fraction of neighbors Qnn(T) which have not changed 
between two images separated by a time interval r, for differ- 
ent packing fractions, as indicated in the legend. 

vious experiments performed with brass discs, we note 
two significant differences. First, the transition is shifted 
to lower values of the packing fraction, a fact which we 
attribute to the higher interparticle friction of the photo- 
elastic discs. Second the range of packing fraction on 
which the pressure varies from say 0.1 to 1 is much larger, 
([0.800 - 0.820] as compared to [0.840 - 0.845]), a direct 
evidence of the discs softness. 

Altogether, the above compression protocol produces 
a structure in terms of nearest neighbor relationships, 
which is completely frozen on experimental timescales. 
The Laguerre tessellation of two different packings sepa- 
rated by a time lag 5000 can be superimposed almost per- 
fectly, even at the loosest packing fraction (fig. [2] bottom- 
left). This is further quantified by Qnn{^)i the average 
fraction of neighbor relationships surviving in a time in- 
terval r. Qnn remains larger than 95% even for the loos- 
est packing fraction and barely departs from 1 for the 
denser ones (fig. [5] bottom-right). 

We now turn to the description of the contact network, 
as provided by the analysis of the photoelastic pictures. 
A threshold is applied to both the gap between neigh- 
boring particles and the contact force, to decide what 
particles are in contact (Fig. [s] (left) ). There is always 
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FIG. 3: Contacts (color online) Left: Map of the coordina- 
tion number for (j) — 0.8125, (grayscale level from black - high 
coordination number - to white - low coordination number. 
Contacts are indicated in red. Right: Average coordination 
number z vs. packing fraction 0, determined according to 
various force thresholds 

some arbitrariness in the choice of such thresholds and 
one should keep in mind that the resulting "contact net- 
work" is actually the union of the real contacts and very 
small gaps. Fig. [3] (right) displays the average number of 
contacts, as a function of the packing fraction for dif- 
ferent thresholds. While the absolute value of z depends 
on the threshold, the dependance on cj) is very robust: z 
is constant at low packing fractions, displays a kink at 
some intermediate packing fraction, and thereafter, it in- 
creases. The kink occurs at a packing fraction which is 
independent of the threshold. We interpret this packing 
fraction as the signature of the crossover to jamming in 
the presence of vibration and find = 0.8143 ± 0.0005. 
Above ^j, the system is jammed and, z ^ a {(j)—(j)j)^-\-Zc^ 
where P G [0.4 - 1] and Zc G [2.7 - 4.3] depend on the 
threshold. Since counting arguments for frictional pack- 
ings [1 constrain Zc between 3 and 4 for two-dimensional 
systems, it is fair to say that the transition indicated by 
the kink is robust to the thresholding procedure. 

The natural quantity of to characterize the dynamics 
of contacts is the contact overlap between t and t -\- r: 



(1) 



where (3f(t,r) = 6(2 - |(5zi(t, r)|), with e(.), the Heavy- 
side function and Szi{t^r)^ the change in number of con- 
tact of grain between t and t -\- r. Fig. [4] (top left) 
displays Qz{t) = {Q^{t^r))t for the various packing frac- 
tions, where {.)t denotes the time average. The black 
curve corresponds to the packing fraction of the jam- 
ming crossover (f)j. For (j) > Q^(j) remains constant 
at values ranging between 0.7 and 0.9, indicating that 
there is no long time decorrelation of the contact net- 
work. The sole decorrelation observed above occurs 
at short time. This decorrelation is more significant when 
(j) approaches and is presumably related to the den- 
sity of rattlers, which also increases. More remarkable 
is the fact that, for (j) < this short times relaxation 
saturates, while a long time relaxation clearly sets in. In 
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FIG. 4: Contact Dynamics (color online). Top Left: Av- 
erage contact overlap Qz{r) for different packing fractions 
(same color code as in fig. [5] bottom right with the curve 
corresponding to plotted in black). Top Right: Con- 
tact overlap dynamical susceptibility xl{^) (same color code). 
Bottom Left: Maximal contact overlap dynamical suscep- 
tibility xV ^) versus and maximal dynamical susceptibil- 
ity X4 C) versus 0, computed from the grain displacements. 
Note the □ data point, which has an anomalously low value: 
we interpret it as the signature of a lack of statistics close 
to the transition, where the timescales of the heterogeneities 
become extremely large. Bottom right: Spatial correlation 
length ^2 of the coordination number versus 0. 

some sense the jamming crossover, signals a transition 
between a frozen and an evolving contact network, the 
phenomenology of which is reminiscent of the glass tran- 
sition scenario [17 , except that it concerns the contact 
network instead of the density profile. 

To further characterize this long term dynamics, we 
compute the dynamical susceptibility (see [18] for an in- 
troduction to dynamical heterogeneities): 



Xl{r)=N 



Var{Q^it,T)) 
{Var{Qf{t,T)))i 



(2) 



where the variances are sampled over time and de- 
notes the average over the grains. This dynamic suscepti- 
bility measures the spatial correlation in the dynamics of 
the contact network. One sees in Fig. |4] (top right) that 
xKt) becomes significant for (j) < and then exhibits 
a maximum xl*? which in turn displays a clear maxi- 
mum at a packing fraction ^* = 0.810 ± 0.001 (Fig. [I] 
bottom left) : the reorganization of the contacts is max- 
imally collective, indicating a dynamical crossover below 
When performing the same experiments with brass 
discs [ill [15], the authors reported a maximum of the 
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correlations of the dynamics, associated with the parti- 
cle displacements. More precisely, they computed x^i^) 
using; formula (H, with Qf(t,r) replaced by Qf(t,r) = 
exp(-||ri(t + r) — ri(t)|p/cr^(r)). Here, fi{t) is the posi- 
tion of grain i at time t, and cr(r) the root mean square 
displacement of the grains. In the present case (Fig. |4] 
bottom left), X4* correlates strongly with xT P^^" 
ticular also has a maximum at Quite remarkably X4* 
is 20 times larger than xl*- This suggest that the correla- 
tion is not trivially a one to one correspondence, and calls 
for further investigations. Furthermore, these dynamical 
heterogeneities find their root in the static heterogeneity 
of the contact network. Indeed, computing the spatial 
correlator of the coordination number (not shown here) 
we extract a correlation length 7 which shares the same 
non-monotonic dependance on (/), with a maximum at (/)* 
(Fig. |4] (bottom right)). 

Let us now summarize and discuss our observations. 
Having prepared a granular glass, the neighborhood 
structure of which is frozen, we have observed two salient 
features of the contact network: (i) the evolution of the 
averaged number of contacts with packing fraction, z{(p)^ 
points to a first transitional packing fraction ^j; (ii) the 
dynamics of the contact network, together with the fluc- 
tuation of the coordination number are maximally het- 
erogeneous at a packing fraction (/)* < 4)j. These results 
clarify those previously obtained in the same setup with 
brass discs ^4]. These authors observed a maximum 
in the heterogeneities of the dynamics for the packing 
fraction, where Pdyn{(I)) and PstatW intersects. They 
were able to attribute their observations to jamming, but 
they could not precisely identify the underlying mecha- 
nism responsible for these heterogeneities. Here, we show 
that they are associated with the dynamics of the contact 
network, which in turn is governed by its heterogeneous 
structure. Note that contacts include both positive over- 
laps and small gaps between particles. 

The existence of a maximal correlation length suggests 
that the experiment is probing both sides of the jam- 
ming transition. In the case of the brass discs, this was 
a puzzling conclusion, given the very large stiffness of 
the discs. Our results resolve this apparent paradox by 
demonstrating that there are several signatures of point 
J at finite vibration, and that the one associated with the 
dynamical heterogeneities occurs at a lower packing frac- 
tion than the one for which the average number of contact 
increases. In fact, in the present study, Pdyn — Pstat 
at 0*. in the case of the brass discs presumably co- 
incides with the very sharp increase of pressure observed 
in [14 at the largest packing fractions. 

Several signatures of point J have also been re- 
ported recently in simulations of thermal frictionless soft 
spheres [6^ 9^ JT . Motivated by the observation that the 
first peak of the partial pair correlation functions exhibits 
a finite maximum at a packing fraction ^j(T) |8j, the au- 
thors of ^ showed that at any finite T (l)j{T) < (l)i{T)^ 



the isostatic packing fraction where z = 2d. Recent dis- 
cussions [19^ confirm that, depending on the observable, 
there are numerous crossover lines meeting at point J 
at T = 0. This phenomenology is reminiscent of the 
so-called Widom lines observed in the supercritical re- 
gion close to a critical point [20l [21] . Our observations 
may find an explanation in this kind of scenario. Alter- 
natively, one cannot rule out the possible effect of static 
friction, which, by adding an extensive number of degrees 
of freedom to the packing, allows for hyper-static pack- 
ings at jamming, and can significantly alter the above 
picture. Future plans are to investigate the role of the 
vibration magnitude to further ellucidate this matter. 
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minating discussions and are grateful to Cecile Wiertel- 
Gasquet and Vincent Padilla for their skillful technical 
assistance. 
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